Rotational spectrum of molecular ion NH"*" as a probe for ol- and me/mp-variation 
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We identify the molecular ion NH+ as a potential candidate for probing variations in the fine 
structure constant a and electron-to-proton mass ratio /i. NH+ has an anomalously low-lying 
excited ''E" state, being only a few hundred cm"^ above the ground '^11 state. Being a light molecule, 
this proximity is such that rotational levels of the respective states are highly intermixed for low 
angular momenta. We find that several low-frequency transitions within the collective rotational 
spectrum experience enhanced sensitivity to ct- and /i-variation. This is attributable to the close 
proximity of the and states, as well as the ensuing strong spin-orbit coupling between 

them. Suggestions that NH^ may exist in interstellar space and recent predictions that trapped-ion 
precision spectroscopy will be adaptable to molecular ions make NH"'' a promising system for future 
astrophysical and laboratory studies of a- and /i-variation. 
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PACS numbers: 06.20.Jr, 06.30.Ft, 33.20.Bx 



I. INTRODUCTION 

The Standard Model of particle physics provides a solid 
foundation from which physical phenomena of strong 
and electroweak nature — from high-energy scattering to 
atomic and molecular structure — can be successfully de- 
scribed. The Standard Model itself does not predict 
precise values of fundamental constants such as the fine 
structure constant a = /Tic or the electron-to-proton 
mass ratio /i = nic/mp, but rather accepts the exper- 
imentally observed values as input parameters to the 
physical theory. Speculative theories which go beyond 
the Standard Model, such as string theories, suggest that 
these constants may vary in time or space [l[, enticing 
both theorists and experimentalists alike to contrive fa- 
vorable means for detecting variations of these constants. 

One method for determining variations — or limitations 
on variations — in a and fi is from analysis of atomic or 
molecular absorption lines originating from interstellar 
space. Comparison with laboratory spectra can, in prin- 
ciple, reveal variations on cosmological time or distance 
scales. Employing this method, groups have reported ev- 
idence for non-zero variations in both a 0, Q and n Q , 
though other analyses have shown no variation at similar 
levels of accuracy [^-Q (see also Comment [l3|)- Comple- 
mentary to the astrophysical studies, terrestrial experi- 
ments probe a- and /^-variation on much smaller time 
and distance scales, but benefit from the high-precision 
and reproducibility that the spectroscopic experiments 
offer. Given in terms of temporal drift, laboratory ex- 
periments have placed stringent, model-free constraints 
on variations of both a [lH and fi [l^ : 



Both the astrophysical and laboratory methods may 
realize significant gains by utilizing atomic or molecu- 
lar species which have enhanced sensitivity to a- and 
/i-variation. Much theoretical effort has been dedicated 
to identifying such systems, with a general strategy be- 
ing to locate accidental near-degeneracies within energy 
spectra. In this spirit, measurements have been p ro p osed 
for multiply-c harg ed ions [l^ [l^l , diatomic (15I - I22I and 
more complex |23l - l26| molecules, and even nuclei |27H3S 
Note that "near-degeneracy" is a relative term here; for 
example, the ground and anomalously low-lying excited 
state of the ^^^Th nucleus — being separated by an inter- 
val of 7.6 eV [13] — are nearly-degenerate relative to the 
typical energy scale of nuclear excitation (> 10 keV). 

The relative sensitivity to a- and /x-variation for a 
given transition may be parameterized in terms of dimen- 
sionless coefficients Qa and Q^, defined by the relation 

Slo ^ ^'^ r\ 

= Qa h Ql3 , 

oj a fj, 

with uj being the transition energy. "Typical" transi- 
tions have sensitivity coefficients on order of unity or less. 
For example, in a diatomic molecule a typical fine struc- 
ture transition has Qa ~ 2, « 0, while a typical 
vibrational transition has Qa ~ 0, w 1/2. Near- 
degeneracies resulting from a cancellation between fine 
structure and vibrational intervals can, however, lead to 
transitions with sensitivity coefficients orders of magni- 
tude larger than unity [Tg. Issj. Several atomic, molec- 
ular, or nuclear transitions which enjoy large enhance- 
ment, however, may be irrelevant for astrophysical stud- 
ies or may prove unfavorable for spectroscopic experi- 
ments, and therefore may be of limited utility. 

In this paper we consider the molecular ion NII+ as a 
candidate for measuring variation in a and fx. NH"'' pos- 
sesses an accidental near-degeneracy between its ground 
with the dot signifying a derivation with respect to time. {X ^H) and first excited (a ^E~) electronic states, these 



a/a = (-1.6 ±2.3) x 10 
fi/n = (3.8 ±5.6) X IQ-^^ 
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TABLE I. Molecular parameters used to describe the rota- 
tional spectrum of the ■^11 and ''E" electronic states of NH+ 
along with their lowest-order scalings with a and fi. Parame- 
ters generally depend on vibrational state v, with some being 
appropriate for both electronic states (e.g., rotational con- 
stants B) whereas others are only appropriate for (e.g., 
spin-orbit constant A) or ^E" (e.g., spin-spin constant A). 
The physical significance of these parameters is described 
in Ref. [iJli the non-standard parameters ^3/21 a-nd 

quantify spin-orbit coupling between ^11 and ''E" states 
of similar « j^,!!!]. 
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being separated by only a few hundred cm~^ (for scale, 
the next electronic state is ~ 22200 cm~^ higher [1^). As 
we will show, this near-degeneracy results in enhanced 
sensitivity coefficients of order 10-100 for a number of 
transitions within the rotational spectrum. Moreover, 
NH+ is a light molecule which has been suggested to be a 
component in interstellar clouds [13, [3^. though to-date 
it has not been detected in such media [S^]- Further still, 
motivated by the search for an electron electric dipolc 
moment, Leanhardt et al. [40| have argued that high pre- 
cision spectroscopy may be performed on molecular ions 
within a Paul trap, with experiments now underway for 
HfF+ [U . These considerations suggest that NH+ could 
serve as a valuable probe of a- and /x- variation in future 
astrophysical or laboratory studies. 

II. ROTATIONAL SPECTRUM OF NH+ 

The rotational spectrum of NH+ was analyzed exper- 
imentally in some detail several years ago by Kawaguchi 
and Amano [l^ and more recently by Hiibers et al. [i^ . 
Both papers tabulate molecular parameters which, to- 
gether with the appropriate effective Hamiltonian, are 
capable of furnishing the rotational spectrum. Evolution 
of the spectrum with respect to a- and //-variation may 
be determined so long as the the scaling of the molec- 
ular parameters with respect to a and // is known. A 
list of molecular parameters relevant to the ^11 and 
states of NH+ are provided in Table U along with their 
lowest-order scalings with a and 

Figure [T] illustrates the rotational spectra of NH+ for 
z; = and v = \ vibrational states. We see that the 
^E^, w ~ state is in such close proximity to the ground 
^n, V — Q state that the corresponding rotational spectra 
begin to overlap for relatively low {N — 5) levels of the 



ladder. For v = 1, the states are closer and even the 
lowest rotational levels are seen to be heavily intermixed. 
This close proximity of and ^Y.^ states allows for siz- 
able coupling between the two via spin-orbit interaction. 
This is accounted for in the effective Hamiltonian through 
the non-standard molecular parameters ^3/21 a-nd 

appearing in Table ID In general, this coupling re- 
sults in noticeable perturbations to the rotational spec- 
trum, with one example being that the ground fi-doublet 
interval is doubled in size j43l |. The effect is more pro- 
nounced for close ^11 and "'S" rotational levels of simi- 
lar angular momentum J and parity p (these being con- 
served quantum numbers). One noteworthy case is the 

2n3/2,u = 0, JP = and ^Y.~,v^ 0,iV = 4, Jp = 
rotational levels (see Fig. [T]). Due to spin-orbit coupling, 
the energy eigenstates are nearly equal admixtures of the 
unperturbed and states. 

In Table |T] we have decomposed the term energy for a 
given electronic and vibrational state into three contri- 
butions, 

T — TiYY ^ ^rcl ^ -^vib 5 

where the respective terms correspond to the non- 
relativistic electronic, relativistic electronic, and vibra- 
tional contributions. For each vibrational subspace, only 
T (^S^) — T (^n) is required to produce the spectrum, 
this being supplied by Refs. [i^, E^. However, to de- 
termine the overall a- and /i-dependence we require the 
partial contributions. The vibrational part may be in- 
ferred from Ref. [l^, and for the w = subspace it is 
found to be 

Tvib - Tvib ('n) = -179 cm-\ (2) 

corresponding to the difference in zero-point energies. 
For the v = 1 subspace this difference is a factor of 3 
larger. The relativistic electronic contribution cannot be 
extracted from experiment; we have determined it by ah 
initio calculation to be 

Trci -Trci fH) =41 cm-i. 

The details of this computation are reserved for the Ap- 
pendix. Constraining r(^S~) — T (^H) to the experi- 
mental value implies a non-relativistic electronic contri- 
bution of 

T„r (^S-) - T„, (2h) = 477 cm-i. 

This last part is insensitive to both a- and /i-variation. 

Finally, in our analysis we choose to neglect the hy- 
perfine splitting of the rotational levels. According to 
Hiibers et al. [43| the largest hypcrfinc constants arc of 
order 100 MHz, or 3 x 10^'^ cm^^. As long as the small- 
est uj considered are of the order 0.1 cm""'^, the hyperfine 
interaction does not change these transition energies by 
more than a few percent. The same must be true for 
sensitivity coefhcicnts Qa and Q^. 
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FIG. 1. Rotational spectrum of NH"*" for (a) d = and (b) v — 1 subspaces. Integer and half-integer labels correspond to 
rotational and total angular momentum quantum numbers and J, respectively, while parity is labeled with ±. Energy is in 
cm^^. 



III. RESULTS 

Before proceeding, we must briefly discuss our man- 
agement of the available experimental data. The ef- 
fective Hamiltonian employed by Hiibers et al. [4^ dif- 
fers subtly from the earlier work of Kawaguchi and 
Amano [42|, with this difference extending to the under- 
lying definitions — and therefore scalings — of the molecu- 
lar parameters found in Table HI For the discussion of the 
previous section we chose to be consistent with Hiibers et 
al., who in turn follow closely the methodical develop- 
ment of the effective Hamiltonian given in the book of 
Brown and Carrington [44| . However, in our analysis we 
have calculated the rotational spectrum and sensitivity 
coefficients by following both references as independently 
as possible, taking care to scale molecular parameters as 
appropriate for each case. In principle, results following 
from Ref. [i^l would be preferred, as the data is more ac- 
curate and the effective Hamiltonian formulation is more 
transparent. Unfortunately, Ref. [13 lacks molecular pa- 
rameters for the V = 1 subspace and in itself does not 
provide information for the vibrational contribution to 
the term energy in the d = subspace [i.e., Eq. In 
this section, we choose to maintain consistency by pre- 
senting our results following exclusively from Ref. [42| . 
Implications of this particular choice will be discussed 
more in Sec. IIVI 

In practice, the dimensionless sensitivity coefficients 



for a given transition arc found from the relation (here 



where uj — Ai? is the energy difference between levels 
and Aqx is the difference between dimensional sensitiv- 
ity coefficients defined by 

«=A-|f. ,4, 

Scaling molecular parameters accordingly, we obtain 
rotational energy levels by diagonalizing the effective 
Hamiltonian for multiple values of fi and a in the neigh- 
borhood of the known prcscnt-day values; numerical dif- 
ferentiation is then used to obtain the coefficients qx for 
each level. Note that proper interpretation of Eq. ^ 
requires us to specify our employed unit system, namely 
atomic units (i.e., the atomic unit of energy 1 Hartree = 
a^rricC? is assumed constant). A more subtle point is 
that the sensitivity coefficients Qx , despite being dimen- 
sionless. also depend on this specification of atomic units. 
Differences Q'x — Q'xi on the other hand, do not depend 
on the choice of unit system. 

Tables [H] - IIVI display our results for select (low- 
frequency) transitions within the NH+ rotational spec- 
trum. Along with the transition energy lo and sensitivity 
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TABLE II. S7-doublet transitions for the v — and v = 1 vibrational states of NH+. The parity p and the energy E are given 
for the upper level of the doublet, with energy being referenced from the lowest level of the respective vibrational subspace. E 
and uj are in cm~^ and ||-E1||^ is in atomic units. 
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coefficients Qa and Q^, we also tabulate reduced matrix 
elements for electric dipole transition amplitudes, 

\\El\\ = i~l) .J, ^ . , 

-M' q M ) 

where Dq is a spherical component of the electric dipole 
operator and M is an angular momentum projection 
along the z-axis in a space-fixed frame. In the non- 
rclativistic limit all ofF-diagonal amplitudes between elec- 
tronic states ^Ili/2, ^113/2, and "'E" vanish. These am- 
plitudes appear only after spin-orbit interaction mixes 
these state. As spin-orbit coupling is embedded in the ef- 
fective Hamiltonian, only the dipole moments of the ^11 
and ^E" states in the molecule-fixed frame are further 
needed to determine the off-diagonal amplitudes. We 
have calculated the dipole moments for these states to 
be 0.91 a.u. and 0.82 a.u., respectively (see Appendix). 



A. r2-doubIet transitions 

Sensitivity of Jl-doublet transitions to a- and li- 
variation has been discussed previously in Refs. |19l - [2l| . 
In NII+ the situation is complicated by the proximity 
and strong spin-orbit coupling of the ^11 and ^E" states, 
resulting in less predictable behavior of the intervals. In 
some instances this leads to anomalously small transi- 
tion energies w, with subsequently magnified sensitivity 
coefficients [see Eq. ([3])]. One example is the J = 15/2 
doublet of the ^113/2, w = state. Here we find that the 
spin-orbit coupling leads to a more than 10-fold decrease 
in the transition energy, with large sensitivity coefficients 
Qa = —82 and = —127 following. Full results for il- 
doublet transitions are given in Table [III 

Another mechanism for enhancement is the state- 
mixing itself, as exemplified by the J = 11/2 doublet 
of the ^113/2, 'y = state. Indeed, spin-orbit coupling has 
a significant effect on this transition energy, even cans- 
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TABLE III. Low-frequency < 30 cm ^) transitions ^11 



Negative frequencies mean that the "^E level is below 
'n level. "E" levels are labeled with quantum numbers 



^E- 
the 

Nj and ^11 levels are labeled with f2j. E and 
and is in atomic units. Transitions with \\E1 
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TABLE IV. Low-frequency {\oj\ < 
^ni/2 ^113/2 ill the V = 1 subspace 
mean that the ^ 



30 cm ) transitions 
Negative frequencies 
Levels are 

labeled with quantum numbers J^. E and uj are in cm~^ and 
WElW^ is in atomic units. Transitions with \\El\f < 10"^ 
a.u. are skipped. 
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ing u! to change sign. This modification of w does not 
promote enhancement, however, as the absolute value of 
Lo is found to increase by about 50%. Enhancement here 
comes rather from the numerator of Eq. ([3]). As men- 
tioned in Secini the negative parity state in this doublet 
is nearly an equal admixture of ^11 and "'E- due to spin- 
obit coupling of near resonant levels. The ^11 portion 
brings no enhancement, while the '^Y." portion provides 
larger Ag^ and Ag^ values expected of ^11 — > tran- 
sitions. In fact, Ag^i may be readily estimated: 



- (-179 cm" 
2 ^ 



1) + (-22) (15 cm-i) 



-210 cm- 



where the first term is from a difference in vibrational 
energies (with a factor of 1/2 from the fi^^^ scaling) and 
the second term is from a difference in rotational energies 
\B w 15 cm-^ for both states and —22 is the difference in 



N{N -1-1)]; the leading 1/2 is a weight factor accounting 
for the fact that the eigenstate is about "half" ^S". The 
comparatively small interval, uj = 3.4 cm-^, implies a 
sensitivity coefficient w —60. This agrees with our 
tabulated result for this transition, wherein Qa = 20 and 
Qa = -53. 



B. 



E transitions. 



From the preceding discussion, we can infer that sen- 
sitivity coefficients Qa and will be enhanced for 
— >■ ^I]- cross transitions having lj comparable to the 
51-doublet intervals (-^1 — 10 cm-^). From the perspec- 
tive of electronic transitions, these are anomalously small 
LJ. There are a number of such transitions in the spec- 
trum, and in Tables Hill we present results for transitions 
with |cj| < 30 cm-^. Again we find sensitivity coefficients 
on the order of 10 — 100. 

Transitions between ^11 and ^I]- states require the 
change of the electronic spin, S* = 1/2 S* = 3/2. Thus, 
we would expect strong suppression of transition ampli- 
tudes. This is found to be generally true for the v ~ 
subspace. The exception is the transition complementary 
to f2-doublet transition described in the previous section 
(involving the other state which is a nearly equal admix- 
ture of and ^S"); this has an "enhanced" transition 
amplitude on the order of the J7-doublet transition am- 
plitudes. For the v — I subspace, the difference in term 
energies is smaller ( — 19 cm"^ compared to 339 cm-^ for 
V = 0) and the mixing caused by spin-orbit interaction 
is generally stronger. Because of this, there are several 
sufficiently strong transitions. 



Fine-structure transitions. 



Fine-structure transitions between II1/2 and TI3/2 
states constitute another source of low-frequency transi- 
tions; following from our analysis of f2-doubling, we may 
suspect enhancement of Qa and here as well. In the 
Hund's case 'a' limit, these transitions require spin-flip 
(i.e.; E = ±1/2 E = T1/2, for spin-projection E on 
the internuclear axis) and are therefore forbidden in the 
non-relativistic limit. The interaction with the "^E" state 
can open some of these transitions. It is a second or- 
der effect and the transition amplitudes are significantly 
smaller than for ^11 ^E- transitions. Table IIVI lists 
a number such transitions within the v = 1 subspace, 
with the most sensitive transition having Qa = 31 and 



-12. 



IV. ACCURACY ASSESSMENT 



As mentioned in Sec. IIIIl we analyzed the NH+ rota- 
tional spectrum following the papers of Kawaguchi and 
Amano [4^ 1 and Hiibers et al. (43j as independently as 



6 



possible. A comparison of these results allows us to as- 
sess the accuracy of our data presented in Tables |lT] - 
IIVI For the majority of the transitions, our values of 
Qa and agree to within a few percent following both 
references. This agreement, however, deteriorates for a 
handful of the transitions, with these transitions tending 
to be the ones with the largest Qa and Q^. The bulk 
of the discrepancy in these cases can be attributed to a 
difference in uj [i.e., the denominator of Eq. ([3|)], which in 
the most extreme case differs by a factor of two. Direct 
measurement of the transition energies can remove this 
uncertainty, provided that the sensitivity coefficients are 
corrected appropriately, 

/ , .thcor. \ 

QT'- = QT°'- [-^) . (5) 

With this prescription, theoretical uncertainty is then 
due to calculated A^q and Ag^. We find that our Ag^ 
are consistent to ^10%, with agreement being somewhat 
better for Ag^ at < 5%. We emphasize that our compar- 
ison here is limited to the w = subspace, as Hiibers et 
al. do not give data for the v = \ subspace. 

The above analysis effectively gauges the spread in Ag^ 
and Ag^ arising from uncertainty in the molecular pa- 
rameters of Table [T] (we note that the difference between 
molecular parameters tabulated in Refs. [13, 13 ^yP" 
ically more than 1-cr, even after the disparity in effec- 
tive Hamiltonians is accounted for H^). The notable 
exclusion is the relativistic electronic contribution to the 
term energy. This is supplied by our computed value, 
Trei {^^~) - Tre\ (^II) = 41 cm'^ for which we estimate 
an uncertainty of 10% (see Appendix). Taking values 
within this window, we find deviations in Aqa to be un- 
der 10% for most transitions, but up to 30% for a select 
few. Since T^^i is insensitive to variations in fi, the Aq^ 
are not affected. 

Finally, we mention that our analysis is based on the 
lowest-order scaling of the molecular parameters. In 
principle, higher-order contributions alter the a- and ^- 
dependence of the parameters. The higher-order contri- 
butions are small, and we expect their effects on Aqa and 
Aq^ to be negligible. 

Altogether, we estimate uncertainty in our Aqa and 
Ag^ to be about 30% and 5%, respectively. With experi- 
mental values of uj and subsequent application of Eq. ([5]) , 
this uncertainty extends to our tabulated sensitivity co- 
efficients Qa and Q^. For most transitions, the 30% es- 
timate is highly conservative. 



V. CONCLUSION 

For laboratory experiments, it is useful to quantify 
absolute shift in addition to relative shift. The largest 
Aqa and Aq^ have magnitudes of about 100 cm~^ and 
400 cm~^, respectively. For temporal variations of a and 
fi at the current laboratory limits, Eq. ([T}, the resulting 



drift in uj is 

\uj\ w (400 cm-i) (4 X 10"" yr'^) « 0.5 Hz/yr. 

Here a-variation is neglected, as it's laboratory con- 
straint is three orders of magnitude tighter than for fi- 
variation. We may conclude that if Hertz-level precision 
is achievable with molecular-ion spectroscopy, then NH"*" 
represents a good system to probe for ^-variation in the 
laboratory. 

The natural linewidth gives a fundamental limit to 
spectroscopic precision. The contribution to the natu- 
ral linewidth from a given decay channel n ^ ti' is pro- 
portional to the product a;'^|(n|D|n')p, with uj being the 
transition energy and (rt|D|n') being the electric dipole 
matrix element connecting the two states. For the v ~ 
subspace, the dominant decay channels are the rotational 
transitions ^ — 1. As the rotational splitting Wrot 
grows with A^, so too does the natural linewidth. For 
large A^, where the Hund's case 'b' limit is realized, we 
find the linewidth to be 

where V is the dipole moment in the molecule-fixed 
frame. With the experimental values of B [i^l and our 
calculated values of V, we establish that Trot ^ 30 Hz for 
rotational levels A^ < 20. For the v = 1 subspace, decay 
to the ground vibrational state opens up as well. The 
vibrational transition energy Wvib is larger than the rota- 
tional transition energies Wrot, whereas the dipole matrix 
element is suppressed in comparison. We find the result- 
ing contribution to the linewidth to be 

X being a factor of order unity which accounts for de- 
pendence of V on internuclear separation R (as with 
T), evaluation of x at the equilibrium separation is im- 
plicit in the first expression). With experimental val- 
ues of Wvib [13 and computed values of x (see Ap- 
pendix), we find Tvib — 20 Hz for the ^H state and 
Fvib — 85 Hz for the state. Assuming practical lim- 
itations (10~'^ — 10~^) X r. depending on statistics, we 
conclude that the natural linewidth should allow mea- 
surement of the transition energies at the required level 
of precision. 

To summarize, we have analyzed the sensitivity of the 
NH+ rotational spectrum to variations in the fine struc- 
ture constant a and electron-to-proton mass ratio ji. We 
find enhanced sensitivity for a number of low-frequency 
transitions within both the w = and v = \ spectra, hav- 
ing sensitivity coefficients on the order of 10 — 100. The 
enhanced sensitivity for these transitions is attributed 
to the near degeneracy of the ground X ^H and excited 
a electronic states and the significant spin-orbit cou- 
pling between them. These results could prove useful in 
future astrophysical and laboratory searches for a- and 
/i-variation. 
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VII. APPENDIX: COMPUTATIONAL DETAILS 

The potential energy curves and the dipole moments 
of the and "^S^ states were obtained using the MOL- 
PRO computational package [45|, within the multiref- 
erence configuration interaction approach with single 
and double excitations (MRCISD). Correlation consis- 
tent Dunnin g au g-cc-pVQZ basis sets were employed for 
both atoms ]4y, |43] ; convergence of calculated spectro- 
scopic constants with respect to the basis set was veri- 
fied. The second order Douglas-Kroll-Hess Hamiltonian 
[isl |49| was used to account for scalar relativistic effects. 
The contribution of spin-orbit coupling was calculated 
by employing the Breit-Pauli (BP) operator. Within 
the MOLPRO package, the lowest-order one- and two- 
electron spin-orbit BP operators are used for comput- 
ing the matrix elements between internal configurations 
(no electrons in external orbitals), while for contributions 
of external configurations a mean-field one-electron Fock 



operator is employed. The error introduced by this ap- 
proximation is generally negligible (soj . 

The calculated potential energy curves were used to ex- 
tract the relativistic electronic contribution to the term 
energy, T^i (^S") -T^ei {^U) = 41 cm'^ Agreement be- 
tween our calculated spin-orbit splitting, 80.4 cm~^, and 
the experimental spin-orbit constant A = 81.7 cm"-'^ [i^ 
gives testament to our accuracy. Our accuracy is likely 
limited by higher-order (^ a^) effects, arising primar- 
ily from the spin-orbit mixing between the ^11 and 
states. This mixing is already explicitly accounted for in 
the effective Hamiltonian by the parameters ^1/2, ^3/21 
and and its effect on the potential energy curves here 
represents somewhat of an intrusion. We estimate these 
effects to be a few cm~^, leading us to ascribe 10% accu- 
racy to our calculated value of Trci {'^'^~) — ^rci (^H) = 
41 cm"^ 

We have computed the dipole moments in the 
molecule-fixed frame (origin at the center-of-mass) to be 
0.91 a.u. and 0.82 a.u. for the ^11 and states, re- 
spectively. Cheng et al. have previously computed the 
dipole moment to be 0.7897 a.u. [5l|, and we suspect 
that their result is more accurate. High-accuracy is not 
required for the dipole moments, and this agreement is 
sufficient for our purposes. We also compute dependence 
of the dipole moments on internuclear separation; for the 
factor X appearing in Eq. we obtain values 0.73 and 
2.0 for -^n and "^E", respectively. 
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